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ABSTRACT 

J  A  modified  form  of  the  Klett  inversion  method  of  the  Ildar  equa¬ 
tion  is  developed.  The  accuracy  of  this  inversion  method  is  explored 
and  validated.  Unlike  previous  inversions,  this  one  does  not  require 
making  a  guess  for  each  lidar  return  to  initiate  the  calculation. 

Included  are  discussions  of  the  other  inversion  methods.  The 
assumptions  contained  in  the  single-scatter  lidar  equation  and  in  the 
Inversion  method  are  outlined  as  is  their  relative  importance.  Special 
attention  is  given  to  the  backscatter  versus  total  extinction  relation, 
multi-scattering  and  non-ideal  detector  response.  _ 


RESUME 

On  a  mod if 16  la  mSthode  d* inversion  d£velopp$e  par  Klett  pour 
1* Equation  lidar.  La  vallditS  et  1' exactitude  de  cette  nouvelle 
m6thode  soot  dSmontrSes.  Contrairoment  aux  mCthodes  d£j£  ccnnues,  la 
modification  perrnet  d'gviter  1* emplol  d'une  estimation  pour  determiner 
1’ extinction  a  cheque  retour  du  rayonnenent  lidar. 

On  discute  aussi  des  autres  mithodes  d' Inversion.  On  soullgne 
1' importance  des  hypotheses  dans  le  diveloppement  del' equation  lidar 
pour  la  diffusion  unique  et  la  m6thode  d' inversion.  On  accorde  une 
attention  sp€ciale  2  la  relation  entre  la  retrodlffusidn  et  1' extinct 
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NOMENCLATURE 


Receiver  area 


B y  Planck  enission  function  at  frequency  v 


C(r)  Clear  air  backscattered  power  at  distance  r 
c  Speed  of  light 

c  As  »  subscript  refers  to  clear  air 

D  Offset  distance  used  in  calculating  digitization 

and  numerical  errors 

d  Constant  of  proportionality  relating  backscatter  and  extinction 

F(r)  System  function  Including  system  sensitivity  and  geometric 
crossover  as  a  function  of  distance  r 


f- 


£ 


H 


Relative  change  between  clear  air  signal  and  peak  of  signal 


Iv  Specific  intensity  of  radiation  of  frequency  v 


Jy  Source  intensity  at  frequency  v 


Power  law  coefficient  relating  backscatter] 


Incremental  distance  through  an  aerosol 


As  a  subscript  refers  to  distance  at  which 


P(r)  Lidar  backscattered  power  return  at  distance  r 


and  extinction 


a  guess  is  made 


f? 

t. 


\ v. 


-W.v.v.^v 
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o  As  a  subscript  refers  to  Initial  distance 

p( y)  Phase  function  as  a  function  of  the  posine  of  the  scattering 
angle  u 

r  Distance  or  range  from  the  Ildar 

S(r)  Logarithmic  range-adjusted  power  at  range  r 

S_  Experimentally  determined  S(r) 

£• 

ST  Theoretically  determined  S(r) 

T  Transmission 

T.  Transmission  as  calculated  by  the  inversion 

inv 

W  Concentration 

z  Parameter  for  adjusting  correction  to  logarithmic  amplifier  and 

multi-scattering 

a  Mass  extinction  coefficient 

0(r)  Backseat ter  extinction  coefficient  at  range  r 
o(r)  Volume  extinction  coefficient  at  range  r 

t  Laser  pulse  duration 

y  Cosine  of  scattering  angle 

v  Frequency 
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1.0  INTRODUCTION 

J 

The  Ildar  equation  Is  the  governing  relation  which  permits  the 
calculation  of  extinction  and  estimates  of  concentration  under  condi¬ 
tions  of  s'ngle  scattering  from  a  given  Ildar  return.  However,  to 
obtain  the  solution,  an  inversion  of  the  equation  must  be  performed. , 

J  s-'qo/-'  ■'«*•*,  >'  j  £***»•■$  . 

£>4c/cSV:l7^>"7  /'  f=ta>Li-  <\Xu.ti„>,  «>•<=*-"/  '  InvarS’**  tf  H* 

Until  Recently,  the  inversion  of  Ildar  returns  to  obtain  extinc-  )  • 

tion  coefficients  or  concentration  profiles  has  been  plagued  by  insta-  I 

bility  and  inaccuracies  (Ref.  1).  The  instability  of  previous 
solutions  is  found  to  be  caused  by  mathematical  and  not,  as  acme  have 
suggested,  physical  reasons.  Inaccuracies  can  occur  by  making  unneces¬ 
sary  assumptions  to  further  simplify  the  equation.  The  solution 
proposed  by  Klett  in  1981  (Ref.  1)  goes  a  long  way  to  improve  this 
situation;  however,  as  will  be  shown  in  this  report,  it  is  not  quite 
ideal.  In  Ref.  2  Lentz  describes  a  way  to  overcame  same  of  these  prob¬ 
lems  at  high  visibilities. 

It  seems  that  the  inversion  of  the  Ildar  equation  (or  radar 
equation  with  attenuation)  has  not  been  glvun  the  attention  it 
deserves.  Indeed,  the  basis  of  the  ideas  used  in  the  recent  solutions 
was  contained  in  a  paper  by  Hltschfeld  and  Borden  in  1954  (Ref.  3). 

This  report  is  an  attempt  to  eliminate  or  significantly  reduce 
the  remaining  problems  while  maintaining  a  simple  algorithm.  With  the 
inversion  of  the  lidar  equation  reduced  to  a  simple,  efficient 
algorithm,  analysis  from  many  lidar  returns  (e.g.  from  the  Laser  Cloud 
Mapper  (LCM),  Ref.  4)  is  possible  in  a  short  time.  Appendix  D  gives  a 
FORTRAN  I V  listing  of  the  final  program  used  in  this  report. 

This  work  was  performed  at  DREV  between  May  1981  and  February 
1983  under  PCN  21B05,  Aerosols. 
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2.0  THE  LIDAR  EQUATION 


To  derive  the  lidar  equation,  It  Is  convenient  to  begin  with  the 
equation  of  radiative  transfer  (Ref.  5): 


dIv 

dir-  -  <rv- V 


[i] 


where  1^  is  the  specific  intensity  of  radiation  of  frequency  v  in  the  r 
direction,  o(r)  is  the  total  extinction  coefficient  and  Jv  is  the 
source  term.  Jv  contains  contributions  to  the  intensity  from  scat¬ 
tering  or  emission  and  can  be  given  by 


Jv“  Bv+  il  P(‘1>  Ivdw  [2] 


where  By  is  the  Planck  emission  function,  p  is  the  phase  function  for 
scattering  and  u  is  the  cosine  of  the  scattering  angle.  Note  that  in 
deriving  [l ]  an  important  assumption  has  already  been  made:  that  the 
distance  between  the  scatterers  is  larger  than  the  wavelength  used.  If 
this  condition  holds,  then  the  scattered  radiation  is  incoherent  and 
noninterfering. 


For  the  next  step  in  deriving  the  lidar  equation,  the  assumption 
is  made  that  Iv  »  Jy.  This  assumes  that  no  significant  amount  of 
emission  or  scattering  from  the  medium  are  placed  in  the  direction  of 
propagation.  Note  that  this  contains  the  critical  assumption  that  only 
single  scattering  events  occur.  In  making  this  assumption  we  obtain 

[3] 


dl 


1  v  -  -o(r) 


Iv  dr 


or 


I..  -  <  a<r')dr' 


V  ov 
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with  I  being  the  Initial  intensity  of  the  beam  before  any  extinction. 
Equation  *  is  of  course  the  well-known  Beer's  law. 


For  light  that  is  being  received  in  the  direction  opposite  tc 

the  beam  direction,  the  intensity  received  will  be  proportional  to  IQV> 

geometric  attenuation  where  A  is  the  receiver  area,  and  the  folded 

rz 

CT 

laser  pulse  length  given  by  — ,  where  c  is  speed  of  light  and  t  is  the 
time  duration  of  the  laser  pulse*  Furthermore,  the  intensity  will 
depend  on  the  probability  of  being  backscattered  at  a  range 
r,  [8(r)  *  p(-l)],  and  finally  a  system  function  F(r)  which  Includes 
factors  such  as  system  sensitivity  and  geometric  crossover  of  the 
receiver  and  laser  beam.  Putting  all  these  factors  together  with  [4] 
and  noting  the  power  received  P(r)  «  I y(r),  we  obtain  the  Ildar 
equation: 


P(T)  -  P„^  * 


-2  lra( r’)dr’ 


where  PQ  is  the  initial  laser  pulse  power.  Here,  the  extinction  term 
in  Beer' 8  law  has  been  squared,  since  the  returned  power  must  traverse 
the  same  medium  twice. 


In  summary,  the  derivation  of  [S ]  has  required  the  following 
assumptions: 


1.  Incoherent  scatterers  (Equation  of  Radiative  Transfer); 


2.  no  significant  emission  Iv  »  Jy; 


3.  first-order  multiple  scattering; 


4.  only  backacatter  is  received,  6(r)  -  p(-l)  [ 0(r)  could  be 
properly  Integrated  over  y}; 
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5.  only  a  plane  wave  la  received,  i.e.  even  light  distribution 
over  area  A  (one  can  use  an  effective  area); 

6.  no  beam  pulse  stretching  (t  does  net  change); 

7.  receiver  encoapasses  laser  beam;  beau  overlap  is  complete 
and  there  is  no  blind  spot  (1/r2); 

8.  invariance  of  the  phase  function  over  small  solid  angles; 

9.  o(r)  is  independent  of  PQ; 

10.  particles  are  not  in  the  shadow  of  each  other  (particles  act 
independently);  and 

11.  the  light  source  is  monochromatic. 

3.0  inversion  OF  THE  LIPAR  EQUATION 

In  the  Ildar  equation  [5],  it  Is  usually  the  term  o(r)  that  is 
the  unknown  of  greatest  Interest.  It  can  be  seen  however  that  this 
equation  has  three  unknowns:  o(r),  0(r)  and  F(r). 

F(r)  is  usually  estimated  or  obtained  by  a  calibration  experi¬ 
ment  reducing  the  problem  to  two  unknowns.  In  order  to  proceed,  a 
further  assumption  must  be  made,  that  is 

S(r)  -  dok  [6] 

Klett  has  made  a  literature  search  which  indicates  that 
0.67  <  k  <  1.0  (Ref.  1).  Indeed,  it  seems  that  for  water  fogs,  k  •  1 
is  the  value  most  often  reported  (with  proportionality  constant 
d  ■  0.05  Sr-1  (Refs.  6,  7,  8)).  However,  for  other  types  of  aerosols 
(e.g.  snow,  burning  phosphorus  clouds,  silica  dust,  etc.),  there  is 
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less  information  available  and  furthermore,  if  the  sire  distribution 
changes  in  regions  involving  the  lidar  signal,  k  will  probably  vary. 

For  a  theoretical  study,  see  Ref.  9.  Mure  will  be  said  on  this  assump¬ 
tion  later  in  this  report  (Chapter  5.0). 

tilth  this  final  assumption  [6],  the  lidar  equation  can  now  be 
solved  for  o(r).  In  order  to  simplify  [5]  the  following  definition  is 
usually  made: 

S(r)  -  ln[r2p(r)]  [7] 


Taking  the  derlvate  of  S(r)  with  respect  to  range, 
lidar  equation  [5]  with  F(r)  a  constant  and  the  assumption 


dS 

dr 


k  do 
o  dr 


-2o 


plus  the 
[6],  we  get 

[8] 


In  deriving  [a],  both  F  and  d  disappear.  This  will  be  important 
in  Section  3.4. 


The  following  three  sections  briefly  review  the  major  inversion 
methods  to  date  and  a  fourth  section  describes  the  development  of  the 
modified  inversion  method.  All  four  sections  discuss  the  strengths  «nH 
weaknesses  of  the  various  methods. 


3.1  Slope  and  Ratio  Method 

This  method  will  only  be  mentioned  briefly  since  its  value  as  a 
solution  has  greatly  diminished  with  the  advent  of  the  never  solu¬ 
tions. 

If  one  assumes  that  over  same  distance  the  aerosol  is  uniform, 
we  can  have  —  ■  0.  Therefore,  from  [8], 
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[9] 

so  that  the  slope  of  a  plot  of  S  versus  r  will  give  us  the  value  of  o 
over  the  homogeneous  region.  This  can  be  a  very  poor  approximation  if 
there  is  noise  or  if  the  cloud  is  inhomogeneous.  In  general,  an 
aerosol  that  produces  a  positive  slope  in  S  (likely  in  a  very 
Inhomogeneous  medium)  will  give  negative  o  by  this  method! 

The  ratio  method  is  just  a  variant  of  the  slope  method  when  it 
is  applied  segment  by  segment  (see  Ref.  10).  Note  also  that  [9]  would 
result  from  the  assumption  that  k  -  0,  i.e.  the  backscatfer  coefficient 
is  a  constant  independent  of  the  cloud  concentration.  Even  with  the 
problems  listed  above,  it  is  the  slope  method  that  baa  been  used  the 
most  often  in  the  past  since  it  is  simple  and  the  only  widely  known 
method. 


3.2  Exact  Forward  Inversion  Method 

In  1954,  Hitschield  and  Bordan  (Ref.  3)  solved  the  differential 
equation  [8]  exactly.  It  is  not  difficult  to  solve  because  it  is  a 
standard  differential  equation,  namely  a  homogeneous  Ricatti  equation. 
The  solution  is 


o(R) 


exp[(S  -  S  )/k] 


(0ol  *  1  /^*P  [<S  ‘  So>/k3dr’> 


[10] 


where  the  subscript  o  refers  to  (he  reference  distance  tq  which  is 
closest  to  the  Ildar  system.  This  solution  requires  knowledge  of 
e*  which  can  be  measured  or  guessed.  However,  as  Klett  emphasized 
(Ref*  1),  this  solution  is  unstable  since  two  relatively  large  terms 
are  subtracted  in  the  denominator  to  obtain  a  small  difference.  Thus, 
if  there  is  an  error  in  do,  the  solution  quickly  becomes  unstable  and 
unrealistic  values  of  c  are  computed  at  modestly  small  r.  See  Refs.  1 


wu  ru  •w^r m*  wu  Vu  «  s*  w sr*  v  **  «-/■**--**  *.-*  su-"  «■  n <**r  *«**  ■**  I’v »  -*  \rm  v.  •. 
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and  2  for  aore  details  on  the  error  analysis  of  this  solution.  It  is 
this  instability  that  lias  caused  the  exact  forward  inversion  method  to 
be  neglected  in  favor  of  the  slope  method. 

3.3  Exact  Backward  Inversion  Method  (Klett^  Method) 


It  seems  surprising  to  find  that,  until  1981,  no  one  had 
published  the  idea  of  making  the  guess  or  measurement  of  0O  in  [lO  ]  at 
the  largest  distance  from  the  lidar  system,  even  though  this  idea  was 
mentioned  in  Ref.  3  and  a  more  generalized  version  was  discussed  in 
Ref.  11.  It  appears  that  same  unsuccessful  attempts  were  made 
(Refs.  12,  13).  However,  using  this  idea,  Klett  obtained  the  following 


solution: 


o(r)  - 


exp  [(S  ~  Sm)/k] 

(of1  +  |  //»  exp  [(5  -  Sm)/k]dr') 


where  the  subscript  m  refers  to  the  reference  distance  r  which  is 

in 

furthest  from  the  lidar  system.  Now  the  two  expressions  are  added  in 

the  denominator,  making  the  solution  more  stable,  and  furthermore,  the 

error  in  a  (which  must  be  estimated  for  each  return)  becomes  less 
m 

important  as  r  decreases  since  the  Integral  in  [ll ]  increases.  Thus 
the  solution  will  tend  to  converge  to  the  correct  o(r)  as  r  gets 
smaller. 


Even  though  this  solution  is  greatly  superior  to  [lO],  some 
problems  still  remain.  If,  for  example,  we  take  the  case  of  high  visi¬ 
bility,  it  can  be  seen  that  the  integral  of  [ll ]  is  very  small  and 
indeed  can  be  considered  near  zero  in  comparison  with  cr  ~ l.  If  this 
happens  we  get: 


o(r)  •  Oa  exp  [(S  -  Sa)/k] 


[12] 
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Thus  all  resulting  extinction  coefficients  becaae  linearly 

related  to  a  and  hence  to  the  error  in  a.  From  this,  it  is  apparent 
a  ib 

that,  in  the  case  of  high  visibility,  the  solution  [ill  is  not  very 
stable. 

This  problem  can  be  overcome  (Ref.  2),  by  calibrating  the  lidar 
system  to  find  the  system  constants.  Indeed,  the  very  need  for  making 
a  guess  (o  )  in  this  solution  is  created  by  the  derivation  of  the  dif- 
ferential  equation  [8]  in  which  the  absolute  constants  drop  out.  The 
guess  can  therefore  be  seen  as  an  attempt  to  put  back  this  information. 
Thus,  in  obtaining  the  system  constant,  one  can  make  a  better  guess. 
After  a  guess  has  been  made,  theoretically  calculated  values  of  S  (S^) 
can  be  computed  from  estimated  extinction  coefficients.  These,  in 
turn,  can  be  compared  with  experimentally  obtained  S  (Sg)  (computed 
directly  from  the  Ildar  return  by  Sg  -  ln(P/C).  The  procedure  in 
Ref.  2  is  to  sum  all  and  match  them  with  Sg  by  making 

Z(SE  -  ST>  "  0  l13I 

The  summation  is  intended  to  minimize  noise  in  the  signal.  This  method 
produces  an  acceptable  inversion  for  the  high  visibility  case. 

3.4  New  Approach  (AGILE) 

A  new  approach  called  AGILE  (Automatically  Generated  Inversion 
of  the  Lidar  Equation)  is  developed  in  this  section.  This  approach  was 
adopted  because  the  following  problems  still  remain  with  the  Klett 
method: 

1.  For  the  high  visibility  case,  several  iterations  are  needed 
to  minimize  the  sus  in  [  13] .  It  is  also  possible  that  this 
condition  is  not  strong  enough.  For  example,  if  one  has  a 
rapidly  varying  cloud  of  low  density,  there  will  be  large 
numerical  errors  occurring  la  the  suimutlon  of  [13  ].  These 
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errors  could  make  the  stmt  zero  without  it  necessarily  being 
a  good  inversion  (see  the  discuss! on  of  numerical  errors  in 
Appendix  C). 

2.  For  the  medlim  visibility  case,  the  convergence  of  the  Klett 
method  is  still  not  fast  enough,  since  perhaps  half  of  the 
values  of  the  lidar  return  are  spent  converging  to  the 
correct  value.  Therefore  only  the  front  half  of  the  shot  is 
good.  A  procedure  of  the  type  used  in  Ref.  2  could  help, 
but  will  require  on  the  average  many  iterations  as  well  as 
the  evaluation  of  the  integral  /*a(r')dr'  with  the  ensuing 
numerical  errors. 

3.  At  low  visibilities,  one  can  expect  rapid  convergence  by  the 
Klett  method,  but  now  another  problem  arises  which  can  make 
the  resulting  Inversion  poor.  This  is  caused  by  the  theo¬ 
retical  limit  of  the  integral  in  [ll  ],  as  will  be  shown 
below.  If  this  limit  is  exceeded  due  to  numerical  or 
experimental  errors  or  incorrect  assumptions,  the  resulting 
inversion  will  give  values  for  a  that  become  increasingly 
smaller  than  the  true  value  as  r  decreases. 


the  proposed  modification  is  an  attempt  to  eliminate  or  signifi¬ 
cantly  reduce  the  above-mentioned  problems  in  the  inversion  of  a  Ildar 
signal. 


If  we  start  with  the  lidar  equation  [5]  and  write  it  for  the 
case  of  the  clear  air  returned  power  value  C(r)  (as  done  in  Ref.  14), 
we  have 


C(r)  -  P0^F(r) 


Add**  “2V 


[*] 


JTJW  V  V\  V  V  SN>  v  .VVN  V.-WtVy  '.T 


r  •-yt.vv  <.  m  .  «r.  v.  r. «r J  «\  o  <1  i/vw  A. 


•t  if 
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where  a  Is  the  clear  air  value  of  the  extinction  coefficient  and  we 
c 

have  tat en  a  constant  over  r.  If  we  now  divide  P(r)  by  C(r)  and 
c 

rearrange,  we  get 


-2a  r 
c 


P(r)  k -2fS(r')dr’ 


This  result  has  assumed  that  the  clear  air  aerosol  and  the  cloud  have 

tne  same  value  of  d  (which  is  not  necessarily  true  but  which  does  not 

affect  the  result).  This  is  because  we  can  take  a  as  an  effective 

c 

clear  air  extinction  coefficient  (since  we  are  not  interested  here  in 
clear  air  extinction).  Taking  the  k-th  root  of  each  side  of  [l5]  and 
integrating  over  r,  one  obtains 


-L  -  1.  o  (fMr’ 
rt  a  rP(r")ik  k  Vr  'r 

A  VcTFT  e 


ro  °<r->« 


"l  r®(r,)dr* 


-  |[l-T(r)2/k]  [16] 


where  T(r)  is  the  transmission  at  the  distance  r.  To  further  simplify 
we  will  assume  that  oc  is  small  enough  to  be  negligible  (i.e. 

Oc  «  k/2r)  to  give 

I  °c  ir<£)1/kdr" "  £l  -  t2/1c1  tl7l 


where  it  is  now  understood  that  P,  C  and  T  depend  implicitly  on  r.  The 
importance  of  this  result  can  be  seen  when  it  is  understood  in  terms  of 
its  physical  significance.  Equation  17  states  that  the  normalized 
Integrated  backscatter  has  a  limit.  In  dividing  the  power  return  P  by 
the  reference  signal  C  (in  this  case,  a  clear  air  shot),  one  calibrates 
the  laser  shot  against  system  constants,  system  noise  and  against  the 
1/r2  factor.  The  resulting  improvement  if  the  signal  is  demonstrated 
in  Appendix  A. 
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The  integral  in  [l7  ]  represents  a  sum  over  distance  or  time  of 

the  power  received  during  one  shot  or  of  the  total  power  received  (or 

equivalently,  the  normalized  integrated  backscatter) .  Equation  17  then 

states  that  this  has  a  maximum,  i.e. 

,  ,  p  1/k 

\  «c  f  <£>  dr-.l.,T.O  [18] 

If  this  limit  is  exceeded  then  there  is  either  significant  experimental 
error  (e.g.  detector  downtime),  numerical/digitization  errors  or  some 
physical  assumption  is  no  longer  valid  (e.g.  multi-scattering).  This 
then  provides  an  Immediate  check,  albeit  weak,  on  the  system,  inversion 
method  and  physical  assumptions.  It  should  be  noted  that  limits  simi¬ 
lar  to  those  in  [l8  ]  can  be  derived  for  larger  valves  of  r  where  is 
no  longer  small  enough. 

Another  feature  of  [l7 ]  is  that,  knowing  the  sua  of  the  power 
received  at  a  given  point,  one  can  immediately  calculate  the  transmis¬ 
sion.  This  allows  knowing  T  without  doing  the  inversion  to  get  o  (and 
consequently  eliminates  another  pass  through  a  numerical  Integration). 


Finally,  the  Integral  in  [l7  ]  is  equivalent  to  the  Integral  in 
Klett's  solution  [ll].  In  fact,  writing  Klett's  solution  [ll  ]  after 
calibrating  with  C(r),  one  obtains 


oc(P/C) 


«c<VC«> 


.  2  1  f  mfP.  ,  , 

+  k  ‘4  <c>  dr 


The  Ildar  equation  can  be  written  as  (from  [l5 ]) 


o  -  <Jc  (P/C) 


V  *• 
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By  rioting  that 


°c<VC.> 


_  2/k 

■  T 


from  [20],  and 


2o  r  .  1/k  2/k  2/k 

-j—  /r"  (§)  dr'  -  (T  “  T.  ) 


a  generalization  of  [ 17 ] ,  it  can  be  seen  that  Klett's  solution  is  just 

a  restatement  of  the  Ildar  equation  (except  that  a  or  T  are 

1&  11 

guessed). 

Again  substituting  [ 17 ]  into  [20 ] ,  we  get  the  solution  in  a 
different  form: 

_ <P./c)1/k _  r,,i 


-1  2  ,r  P 

ac  -  4  /:  <!>  ** 


This  is  similar  to  the  forward  inversion  method  (Section  3.2)  except 
that  0Q  is  replaced  by  0£.  Or,  equivalently,  Klett’s  solution  becomes 


(P/C)1'  * 
j.  2  (X  *P* 


/«.  +  t  /;*%>  dr' 


Equations  23  and  24  are  equivalent  since  they  guarantee  that 


STCO  •  Sg(r) 


at  each  Interval  of  r,  a  condition  much  stronger  than  [l3],  providing 
that  errors  do  not  Invalidate  [ 18  ]  •  This  can  be  seen  because  of  [l7] 
and  [20]  and  the  definition  of  and  Sg.  This  condition,  [25],  is 
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noC  satisfied  by  Klett's  solution  [ll]  because  of  possible  errors  in 

the  guess,  a^.  With  solutions  [23]  or  w.  once  a  calibration  shot 

has  been  made  and  a  value  of  a obtained,  no  guessing  will  be  needed  as 

long  as  the  value  of  a^  and  the  calibration  remain  valid.  (The  values 

of  o  can  be  obtained  by  a  trial  inversion  of  one  shot  in  which  the 
c 

transmission  is  measured,  or  by  other  means).  This  has  great  advan¬ 
tages  over  the  other  methods  when  many  shots  under  the  same  conditions 
must  be  processed.  Plus,  the  necessary  condition  [25]  is  automatically 
satisfied,  whereas  in  the  Klett  method  the  weaker  condition  [13]  must 
be  converged  to  and  may  still  give  significant  errors. 

4.0  VALIDATION  OF  AGILE 

In  this  chapter,  evidence  of  the  success  of  AGILE  will  be 
reviewed  and  compared  with  Klett's  solution.  The  remaining  diffi¬ 
culties  will  be  discussed,  as  will  an  Indication  of  the  numerical  and 
digitization  errors.  The  Ildar  data  used  In  this  section  was  obtained 
from  the  DREV  Laser  Cloud  Mapper  (LCM)  (Ref.  4),  with  clouds  produced 
by  burning  :d  phosphorus. 

What  constitutes  a  good  Inversion  of  the  Ildar  equation  [5] 
along  with  the  assumption  [6]? 

It  Is  necessary  and  sufficient  that  a  good  inversion  satisfy  the 
condition  [25]  and  uses  the  correct  value  of  0c*  Equation  13  Is  only  a 
necessary  condition. 

Since  many  assumptions  have  been  made  to  obtain  [23]  or  [24], 
how  can  a  check  be  done  to  indicate  that  the  results  from  the  inversion 
are  a  good  representation  of  the  real  values?  The  following  items  will 
be  compared  and  discussed  in  detail  in  this  chapter  in  order  to  indi¬ 
cate  the  successes  and  limitations  of  the  current  inversion: 
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1.  the  transmission  of  this  inversion  with  respect  to  that 
obtained  from  Klett's  method  (with  the  estimation  procedure 
as  outlined  in  Ref.  2); 

2.  the  transmission  with  that  obtained  directly  by  the  LCM; 

3*  the  extinction  coefficients  and  concentrations  found  typical 
by  this  method  in  burning  red  phosphorus  smoke  and  those 
from  in  situ  measurements; 

4.  cloud  extent; 

5.  clear  air  extinction  coefficients  in  front  of  and  behind  the 
cloud;  and 

6*  the  results  obtained  should  be  Independent  of  the  inversion 
starting  points  in  the  shot. 


4.1 


Vlth  this  new  approach,  two  calibration  shots  are  needed:  one 

in  clear  air  (or  through  a  uniform  medium)  and  the  other  with  same 

transmission  loss  in  order  to  determine  a  .  With  the  LCM  this  is  not  a 

c 

problem,  since  this  system  is  capable  of  making  100  shots  per  second, 
all  in  different  directions  and  many  through  clear  air. 


Once  and  a  good  dear  air  shot  have  been  obtained,  the  big¬ 
gest  remaining  problem  results  from  the  possible  errors  in  the  evalua¬ 
tion  of  the  Integral  in  [23]  or  the  Integral  and  Ta  in  [24].  This  is 
caused  mainly  by  errors  arising  from  the  system  response  or  multi- 
scattering.  For  if  the  system  response  is  too  slow  in  recovery  and/or 
the  cloud  is  dense  enough  for  multi-scatterlcg  to  be  present,  the  limit 
[l8]  could  easily  be  surpassed.  If  this  happens,  neither  [23]  nor  [24] 


•• .  •  \  .  •  vw-y  •  v.y .  ••  v-  v.\  ■  vv. 
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can  work,  and  neither  can  the  Klett  solution.  This  occurs  simply 
because  the  lidar  equation  has  become  invalid. 

If,  however,  the  system  response  is  not  too  bad,  a  correction 
can  be  applied.  Multi-scattering  poses  a  much  more  serious  problem  and 
there  are  currently  no  simple  or  accurate  corrections  for  this  phenome¬ 
non.  Therefore  it  is  suggested  that  in  performing  an  inversion  using 
the  lidar  equation  [5],  a  guarantee  of  the  limit  [l8  ]  should  be  given. 
This  can  be  done  by  decreasing  the  return  values  in  the  individual 
shots  as  the  range  increases,  according  to  how  rapidly  the  normalized 
integrated  backscatter  (or  limit  [l8])  is  approached.  Since  the  func¬ 
tional  form  needed  for  reducing  the  return  values  is  not  generally 
known,  estimates  must  be  used  (which  adds  some  uncertainty  to  the 
results).  However,  one,  needs  to  know  this  function  only  in  the  extreme 
case  of  very  dense  clouds  and  can  be  calibrated  to  some  extent. 

The  functional  form  used  for  the  multi-scattering  is  the  same  as 
that  used  to  correct  for  the  logarithmic  amplifier  (see  Appendix  B). 
Thus  a, single  correction  term  (or  the  same  form  used  by  Ref.  15)  with 
one  adjustable  parameter  attempts  to  correct  for  multi-scattering  and 
the  amplifier  response  (and  any  other  reason  that  may  cause  the  limit 
in  [l8]  to  be  exceeded).  This  factor  is  just  T  *  with  z  being  the 
sd justable  parameter.  Since  both  effects  are  combined,  z  will  be 
larger  than  if  only  one  effect  were  present. 


With  the  data  presented  in  this  report,  z  ■  0.8  was  found  to 
give  good  agreement  with  target  reflectance.  The  value  of  this  para¬ 
meter  in  Ref.  15  is  shown  to  be  0.2  <  z  *  0.3,  with  pure  multi- 
scattering,  a  field  of  view  between  5  and  10  mrad  and  a  size  dlstribur 
tlon  similar  to  that  considered  in  the  above-mentioned  data.  Thus  for 


the  value  of  s  -  0.8,  about  30Z  is  caused  by  multi-scattering.  This 
factor  is  used  only  when  the  optical  depth  becomes  less  than  unity. 


SENT 
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FIGURE  3  -  Comparison  of  transmissions  by  Inversion  versus  by  targets 


4.2  Transmission  Comparisons 

A  comparison  between  the  transmission  computed  from  the  present 
method  and  that  computed  by  Klett's  method  (Ref.  1)  Is  shorn  In  Fig.  1. 
The  graph  shows  that  there  la  excellent  agreement  between  the  two 
methods.  However,  there  Is  still  soms  scatter.  One  reason  for  dis¬ 
crepancies  Is  the  difficulty  In  making  the  proper  guess  In  the  klett 
method  In  denser  clouds,  especially  when  there  is  significant  signal  at 
the  point  In  the  shot  where  the  guess  Is  being  made.  This  difficulty 
Is  due  to  numerical  and  digitisation  errors  that  occur  when  calculating 
the  Integrals  and  obtaining  the  transmission. 
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Another  source  of  error  is  the  response  of  the  amplifier  used  in 
processing  the  return  signal.  Figure  2  shows  the  response  of  the 
logarithmic  amplifier  in  the  LCM.  The  vertical  scale  is  the  loga- 
rlthsLic  output  in  arbitrary  units,  and  the  horizontal  scale  is  the  time 
in  nanoseconds.  Froa  this,  and  with  a  digitization  rate  of  10  ns,  one 
can  see  that  the  rise  tine  of  this  amplifier  is  adequate  hut  that  the 
downtime  will  cause  problems ,  especially  when  the  signal  is  large  and 
varies  rapidly.  Thus,  in  attempting  to  correct  for  this  response  (see 
Appendix  B  or  Section  4.1),  errors  will  occur  in  both  programs  since 
perfect  correction  is  impossible. 

Figure  3  is  a  plot  of  the  transmission  calculated  from  the  new 
algorithm  versus  the  transmission  as  measured  by  the  intensity  of  the 
signal  returned  by  cooperative  targets.  Again  the  agreement  is  reason¬ 
able  but  there  is  same  discrepancy.  Most  of  this  discrepancy  can  be 
explained  by  the  variation  in  target  uniformity  and  errors  in  the  laser 
firing  angle.  These  errors  account  for  a  variation  of  about  ±  20Z. 

Other  errors  that  can  cause  disagreasmnt  in  Fig.  3  are  the 
problems  outlined  above,  i.e.  amplifier  response,  numerical  and  digiti¬ 
sation  errors,  plus,  perhaps,  significant  multi-scattering  for  trans¬ 
mission  such  less  than  402  (about  one  optical  depth).  A  least  mean 
squares  fit  of  the  points  gives  a  straight  line: 

Tinv  "  °’91TE  +  f26l 

with  a  correlation  coefficient  of  0.81,  where  T^Qy  and  Tg  are  the 
tranamiselooe  obtained  from  inversion  and  experiment,  respectively.  The 
correlation  coefficient  value  seems  most  reasonable  when  tt» . variation 
in  the  target  reflectance  la  considered. 
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4.3  Concentration 

As  already  mentioned,  the  data  here  were  taken  by  the  LCM  from  a 
red  phosphorus  generated  smoke.  By  using  values  for  the  mass 
extinction  coefficients  (taken  from  the  DREV  silo  and  9  50Z  relative 
humidity) ,  a  conversion  of  the  total  extinction  coefficient  into  con¬ 
centration  is  possible  by  the  relation 

a  -  oW  ,  [27] 

where  a  is  the  mass  extinction  coefficient,  V  is  the  concentration  in 
grams  per  cubic  centimeter  and  L  is  the  path  length  in  centimeters. 

This  relationship  has  been  verified  by  many  studies  (both  theoretical 
and  experimental)  over  widely  varying  types  of  aerosols,  densities  and 
wavelengths  (see  for  example  Refs.  16-19).  With  such  a  calculation  of 
concentration,  a  comparison  between  the  calculated  and  measured  concen¬ 
trations  can  be  made.  Unfortunately  no  direct  measurement  of  the  con¬ 
centration  was  made  when  these  Ildar  data  were  taken,  but  the  purpose 
here  is  to  show  that  the  concentration  values  obtained  by  the  inversion 
are  reasonable  '‘ballpark’’  figures. 

DREV  and  others  (Refs.  16  ,  20-22)  indicate  a-  1.5-2. 5  m2/g  at 
1.06  pm  and  50Z  R.H.  Furthermore,  one  can  obtain  from  the  literature 
In  situ  field  values  of  the  concentration  for  red  phosphorus  generated 
smokes  (Ref.  23).  Table  I  gives  a  comparison  between  typical  varia¬ 
tions  in  concentrations  between  the  calculation  In  this  report  and , in 
Ref.  23.  The  correspondence  between  the  two  indicates  a  satisfactory 
agreement.  Thu*  the  calculations  dona  here  give  a  good  order  of  magni¬ 
tude  estimate  of  the  actual  field  values  during  the  experiment. 


CLOUD  WIDTH  BY  INVERSION  (M) 
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TABLE  I 


Coaparison  of  field  trial  concentrations  of 
_ red  phosphorus  saoke 


ag/a3 

Smoke  Week  III 

5-200 

(Ref.  6) 

Present  Study 

0.1-300 

CLOUO  WIDTH  (M) 

FIGURE  4  -  Coaparison  of  cloud  extent  calculated  by  Inversion  and  by 
photographic  aeasureaent 


Example  of  raw  Ildar 
RP  smoke 
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4.4  Cloud  Extent 


During  the  trials,  photographs  were  taken  of  both  the  front  and 
back  of  the  cloud.  This  allows  one  to  estimate  the  cloud  extent  photo¬ 
graphically.  Figure  4  is  a  plot  of  the  calculated  cloud  extent  along 
one  laser  line  and  the  estimate  from  the  photographic  data.  Differ¬ 
ences  can  obviously  occur  because  of  the  difficulty  in  estimating  the 
location  of  the  cloud  edges.  In  the  photograph,  this  is  due  to  lack  of 
contrast  and  in  the  lidar  data,  it  is  due  to  the  definition  of  what  is 
cloud  and  what  is  clear  air.  The  agreement  is  good.  The  importance  of 
this  agreement  will  be  discussed  in  the  next  chapter  in  relation  to  the 
assumptions  made. 

4.5  Clear  Air  Extinction 

Figures  5  and  6  are  examples  of  a  raw  lidar  signal  and  fie 
resulting  inversion  for  one  shot.  Figure  6  shows  the  calculated 
extinction  coefficient  plotted  against  distance  from  the  LCM  (or  lidar 
system).  The  important  feature  to  notice  here  is  that  the  clear  air 
extinction  values  in  front  of  and  behind  the  cloud  agree  very  well. 

This  is  essential  and  verifies  the  calculated  transmission  through  the 
cloud. 


4.6  Invariance  of  Inversion  Point 

A  necessary  property  that  all  inversions  should  have  is  that  the 
inversion  should  be  Independent  of  the  initial  starting  point,  l.e. 
overlapping  results  of  differently  initiated  inversions  should  be  iden¬ 
tical  .  The  new  solution  guarantees  this  because  of  the  condition 
Sg  •  However  for  the  Klett  method,  errors  in  the  guessed  value 
will  in  general  produce  values  different  from  guesses  made  at  other 
distances.  Thus  the  resulting  calculation  will  be  different.  As  the 
cloud  density  increases  and/or  varies  rapidly,  this  problem  becomes 
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nore  Important,  especially  for  guesses  made  Inside  the  cloud.  This 
invariance  is  obtained  because  of  the  normalization  of  the  lidar  return 
by  the  clear  air. 

5.0  QUALITY  OF  ASSUMPTIONS 

Many  assumptions  have  been  made  in  the  derivation  of  [23]  and 
[24]  and  their  subsequent  use.  An  attempt  is  made  here  to  identify  the 
most  important  assumptions  and  their  effect  on  the  accuracy  of  the 
inversion.  Most  of  the  assumptions  leading  to  the  lidar  equation 
(Chapter  2.0)  are  excellent  approximations  to  the  conditions  under 
which  the  current  data  were  taken.  Major  exceptions  are  the  assumption 
of  single  scattering  (for  dense  clouds)  and  the  invariance  of  the  phase 
function.  In  this  case,  it  is  a  possible  variation  in  the  size  dis¬ 
tributions  of  the  particles.  This  latter  assumption  is  also  involved 
in  the  relation  [6]  and  when  k  •  1,  which  is  a  constant  with  respect  to 
r.  But  perhaps  Che  most  important  effect  is  that  of  the  logarithmic 
amplifier's  response. 

5.1  Amplifier  Response 

The  time  response  of  the  amplifier  has  already  been  given  in 
Fig.  2.  A  good  check  to  see  if  the  response  of  the  mpllfler  is  being 
adjusted  for  correctly  (see  Appendix  B)  is  to  pass  a  lidar  return  from 
a  target  through  the  inversion  with  this  correction.  Figure  7  is  the 
result  of  such  a  calculation.  Note  that  the  program  has  now  considera¬ 
bly  compensated  for  the  downtime  and  a  relatively  sharp  response 
remains. 

The  target  response  can  be  considerably  enhanced  with  a  Fourier 
transform  and  deconvoluted  with  a  Wiener  filter  (Ref.  24),  but  the 
results,  as  Indicated  In  Fig.  8,  show  that  the  noise  is  considerably 
amplified  and  the  absolute  values  become  unreliable. 
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Figure  9  is  an  illustration  of  results  obtained  from  a  scan  of 
an  HC  (hexachloroe thane)  smoke  without  compensation  for  the  logarithmic 
amplifier.  The  disturbance  in  this  cloud  map  occurs  for  extinction 
coefficients  as  large  as  10“  3  m-1.  Figure  10  shows  the  same  cloud  once 
the  logarithmic  amplifier  response  has  been  corrected.  Notice  that 
even  at  the  lowest  extinction  coefficient,  which  is  considerably 
smaller  than  any  in  Fig.  9,  there  are  no  distortions  of  the  cloud  con¬ 
touring.  The  cloud  maps  of  Figs.  9  and  10  axe  computer  calculated  and 
drawn  and  are  in  the  form  seen  by  the  user  of  the  program  (Ref.  25). 
Note  that  similar  distortions  due  to  multi-scattering  could  occur  in 
clouds. 
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5.2  Multi-Scattering 

The  various  experiments  and  theory  that  have  so  far  been  con¬ 
sidered  in  the  literature  make  it  very  difficult  to  ascertain  the 
effects  of  multi-scattering  for  a  given  beam  geometry,  field  of  view, 
Ildar  system  distance,  aerosol,  wavelength,  size  distribution,  etc. 
(Refs.  15,  26-28).  Indeed,  many  of  the  parameters  on  which  multi- 
scattering  depend  may  only  be  known  vaguely  or  not  at  all  (e.g.  parti¬ 
cle  size  distribution  and  particle  shape). 

When  the  cloud  is  dense  enough,  significant  multi-scattering 
will  invalidate  the  Ildar  equation.  Thus,  either  a  new  Ildar  equation 
Incorporating  multi-scattering  should  be  used  or  the  data  adjusted  so 
that  only  the  single-scatter  component  of  the  signal  is  considered. 

By  examination  of  the  available  theories  and  experiments  (e.g. 
Refs.  15,  26-28),  a  good  rule  of  thumb  would  be  to  consider  multi- 
scattering  when  the  optical  depth  is  greater  than  0.1  for  Ildar  systems 
approximately  1  km  or  more  from  the  cloud  and  perhaps  for  optical 
depths  greater  than  1.0  when  the  Ildar  system  distance  is  dose  (about 
50-100  m),  as  is  the  case  with  the  LCM.  This  assumes  that  the  field  of 
view  of  the  Ildar  system  is  roughly  in  the  range  of  1  to  10  mrad. 
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The  program  that  produced  the  examples  for  this  report  using  the 
new  Inversion  method  guarantees  that  the  Integral  limit  [ 18 ]  will  not 
be  exceeded  so  that  some  compensation  for  multi-scattering  will  be 
present.  No  claim  la  made  here  that  this  compensation  is  complete  or 
correct. 

5.3  Particle  Sire  Distribution 


It  Is  difficult  to  Imagine  conditions  in  nonsolid  particulate 
clouds  or  smokes  in  which  the  particle  size  distribution  does  not 
change.  Varying  humidities  will  affect  the  size  distribution  of  hygro¬ 
scopic  particles  for  example.  For  burning  red  phosphorus,  the  size 
distribution  will  become  stable  only  at  some  distance  from  the  source, 
even  if  the  tumidity  is  uniform.  However,  it  was  shown  in  Section  4.4 
that  there  is  good  agreement  between  the  calculation  and  the  photo¬ 
graphic  evidence.  Thus  the  value  of  k  used  cannot  vary  appreciably 
from  unity  since  this  would  change  the  cloud  extent  and  hence  bias 
Fig.  4,  Since  the  inversion  would  give  larger  values.  This  can  be 
clearly  seen  in  Figs.  11a  and  b,  which  show  a  plot  of  extinction 
coefficient  versus  distance  for  various  values  of  k.  Values  of  k  less 
than  unity  affect  the  cloud  extent  to  a  greater  degree  than  k  >  1. 
Notice  that  the  extinction  values  are  also  changed,  which  would  affect 
the  relationships  already  shown  in  Fig.  3  and  Table  I. 

A  survey  of  the  literature  indicates  that  for  most  conditions, 
k  -  1.  Theoretical  results  (Refs.  7,  29-32)  and  experimental  results 
(Refs.  8,  17,  33-35)  over  a  variety  of  materials,  particle  shapes,  and 
extinction  coefficients  indicate  that  this  assumption  of  linearity  is 
excellent.  However  caution  is  required  since,  as  indicated  by  results 
in  Refs.  8,  9,  33,  36-37,  effects  such  as  multi-scattering,  changes  of 
size  distribution  (due  to  humidity  or  settling,  etc.)  and  index  of 
refraction,  and  particle  shape  can  cause  nonlinearity.  Values  of 
k  -  0.66  have  been  reported  (see  Ref.  31  and  the  references  contained 
therein) . 


DISTANCE  CIO 


FIGURE  11a  -  Effect  of  power  law  variation  In  the  backscatter- 
extinction  relation  on  Inversion  results  (k  <  1) 
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As  mentioned  in  Section  3.4,  the  clear  air  extinction  coeffi¬ 
cient  a  can  be  Interpreted  as  the  effective  extinction  coefficient  in 
order  to  account  for  the  differences  between  the  proportionality  con¬ 
stant  d  in  [6]  relating  the  clear  air  aerosol  to  the  aerosol  under 
study.  The  subsequent  Inversion  using  the  calibration  value  of  oc  will 
continue  to  be  valid  (under  conditions  discussed  in  Section  4.1)  unless 
the  constant  d  in  [6]  for  the  aaterlal  studied  changes. 


6.0  SUMMARY  AND  CONCLUSIONS 


In  this  report,  a  method  has  been  developed  that  permits  the 
Inversion  of  the  Ildar  equation  under  a  vide  variety  of  conditions. 

This  solution  overcomes  instability  in  both  the  high  and  low  visibility 
cases,  in  contrast  with  the  solutions  proposed  by  Klett  and  Lentz 
(Refs.  1,  2),  by  normalizing  the  signal  by  a  clear  air  return  and 
considering  multi-scattering.  This  solution  is  designed  to  invert  the 
Ildar  equation  norm  accurately  and  with  less  computation. 


It  is  demonstrated  that  there  is  a  theoretical  limit  to  the 
total  power  received  (normalized  Integrated  backscatter)  in  a  return 
signal  which  provides: 

1.  a  check  on  the  validity  of  the  Ildar  equation; 

2.  a  check  on  the  system  Itself; 

3.  transmission  at  any  point.  Independent  of  the  inversion; 

4.  that  the  transmission  is  obtained  after  only  one  numerical 
integration  instead  of- two  and  is  therefore  more  accurate, 
as  a  consequence  of  #3;  and 
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5.  a  new  inversion  nethod  without  the  need  for  a  guess  or  a 
need  for  iteration  to  yield  convergence  to  soae  constraint 

The  solution  is  validated  in  a  nuaber  of  ways,  including: 

1.  comparison  with  Klett's  aethod  (Refs*  1  and  2); 

2.  comparison  of  predicted  and  aeasured  transaissions; 

3.  coaparison  of  calculated  concentrations  with  aeasured  con¬ 
centration  values  f roa  the  literature;  and 

4.  cloud  extent. 

It  is  also  shown  that  calibrating  the  Ildar  signal  frith  a 
reference  signal  (in  this  report,  a  clear  air  signal)  aakes  the 
r2  correction  unnecessary  and  considerably  iaprovea  the  ability  to 
detect  a  weak  signal  aaongst  significant  systea  noise.  As  a  result, 
the  values  of  extinction  produced  by  the  inversion  are  independent  of 
the  starting  point  of  the  inversion. 
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APPENDIX  A 

The  Effects  of  Clear  Air  Calibration 

In  Section  3.4,  where  the  nev  inversion  method  is  developed,  the 
Ildar  return  signal  is  divided  by  a  clear  air  return  (or  by  any  other 
reference  signal).  In  so  doing,  the  data  are  immediately  corrected  for 
1/r2  attenuation,  the  system  constant  F(r),  and  most  system  noise.  In 
this  appendix,  the  sensitivity  gained  due  to  the  latter  effect  will  be 
discussed . 

Two  histograms  are  given  in  Figs.  A1  and  A2  of  the  same  Ildar 
return  showing  the  distribution,  after  Inversion,  of  extinction  coeffi¬ 
cients.  Here  Fig.  A1  la  the  Inversion  w‘vh  just  the  1/r2  correction, 
while  Fig.  A2  Is  the  same  Inversion  but  using  the  clear  air  return  for 
calibration.  In  both  cases,  the  shot  has  been  corrected  for  laser 
power  signal  with  respect  to  the  calibrating  shot.  Also,  the  clear  air 
extinction  coefficient  was  taken  to  be  2.0  x  10”5  m~l.  For  this  case, 
the  standard  deviation  is  essentially  halved  by  using  the  calibrating 
shot.  And,  significantly.  Fig.  A2  allows  for  the  detection  of  a  very 
weak  return  of  the  cloud  while,  from  the  statistics  of  the  noise. 

Fig.  A1  does  not.  Indeed,  the  one  point  In  Fig.  A2  (at  2.6  x  I0”5  m”1) 
shows  that  the  cloud  Is  more  than  eight  standard  deviations  from  the 
mean  whereas  it  Is  only  three  standard  deviations  from  the  mean  in 
Fig.  A1  (at  2.3  x  10~5  m“l). 

Figure  A3  la  the  distribution  for  the  extinction  values  of  over 
1200  clear  air  extinction  coefficients.  In  this  way,  an  excellent 
value  of  the  standard  deviation  of  the  clear  air  extinction  In  the  data 
can  be  derived,  enabling  the  detection  of  very  weak  returns.  For 
example.  If  an  extinction  value  Is  over  three  standard  deviations 
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FIGURE  A3  -  Distribution  of  over  1200  extinction  coefficients  with 
deer  eir  normalization 
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APPENDIX  B 


Correction  for  the  Logarithmic  Amplifier  Response 

As  shown  in  Fig.  2,  the  logarithmic  amplifier  response  is  far 
from  ideal,  and  the  use  of  digital  filters  can  be  difficult  (as  in  the 
one  attempt  discussed  in  Section  5.1).  The  method  used  in  the  current 
inversion  program  to  correct  for  this  response  is  to  multiply  the 
return  signal  by  a  factor  less  than  unity  if  the  cloud  is  sufficiently 
dense.  This  factor  decreases  continuously  a8  the  integrated  back- 
scatter  gets  larger.  The  actual  factor  used  is  T  z  where  T  is  the 
optical  depth  and  z  is  the  adjustable  parameter. 

The  program  halts  if  too  much  correction  is  required,  assuming 
that  calculation  beyond  this  point  produces  meaningless  information. 
The  amount  of  correction  needed  and  when  to  halt  the  program  must  be 
empirically  decided.  This  is  done  by  comparison  with  transmission 
values  and  cloud  extent.  Note  that  the  amount  adjusted  for  is  con¬ 
trolled  by  a  single  parameter  z  and  that  this  is  made  a  constant. 

A  multi-parameter  model  was  felt  to  be  too  flexible,  i.e.  it  could  be 
made  to  fit  almost  any  condition.  For  the  program  used  here,  a  value 
of  z  ■  0.8  waa  found  to  be  reasonable;  that  is,  it  was  a  good  fit  with 
the  observed  transmissions. 
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APPENDIX  C 


Numerical  and  Digitization  Errors 


In  calculating  the  inversion,  a  numerical  integration  oust  be 
performed,  with  its  inherent  errors.  Furthermore,  the  digitization  of 
the  analog  signal  obtained  by  the  LCM  also  places  an  unavoidable  error 
in  the  data,  which  is  due  to  the  finite  word  size  of  the  digitizer. 

The  purpose  of  this  appendix  is  to  estimate  the  importance  of  these  t 
errors.  The  errors  produced  in  the  Integration  and  digitization  of  t..e 
artificial  signal  shown  in  Fig.  Cl  will  be  discussed.  This  is  mathe¬ 
matically  a  simple  representation  of  a  range-corrected  Ildar  return 
sijgnal,  with  the  height  H  and  an  offset  D  being  adjustable.  The  para¬ 
meter  H  is  taken  as  a  relative  change  between  the  clear  air,  signal  and 
the  plateau  height  and  D  is  in  units  of  the  digitization  rate. 


The  error  will  be  taken  here  as  the  difference  between  the  theo¬ 
retical  area  under  the  curve  (to  the  first  data  bln  after  the  signal 
height  H)  and  the  numerically  derived  area. 


SIGNAL 


(DISPLACEMENT) 


FIGURE  Cl  -  Artificial  Ildar  return  used  for  estimating  numerical  and 
digitisation  errors 
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TABLE  Cl 

Variation  of  Integration  Error  with  Offset  and 
Signal  Height  (Infinite  Word  Site) 


Signal  Height  (H) 

1 

10 

100 

1000 

0 

0 

<0.5 

<0.5 

<0.5 

0.1 

0 

3 

8 

10 

0.5 

0 

13 

55 

90 

0.9 

0 

30 

150 

470 

OFFSET  (D) 


TABLE  CII 

Variation  of  Digitisation  Error  with  Offset  and 
_ Signal  Height  (8-Blt  Word  Size) 


Signal  Height  (H) 

1 

10 

100 

1000 

0 

0 

1 

50 

90 

OFFSET  (D) 

0.1 

0 

20 

80 

110 

0.5 

0 

30 

120 

300 

0.9 

0 

40 

170 

800 

TABLE  CIII 

Variation  of  Digitisation  Error  with  Offset  aicd 
_ Signal  Height  (16«Blt  Word  Slse) 


OFFSET  (D) 


Signal  Height  (H) 

1 

10 

100 

1000 

0 

0 

0.8 

8 

10 

0.1 

0 

13 

40 

70 

0.5 

0 

20 

90 

200 

0.9 

0 

30 

ISO 

\v'.v 


.*.Yv  .• 
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Cl . 1  Integration  Errors 

Both  the  offset  D  end  height  H  of  e  'Idar  return  can  vary  and. 

In  the  case  of  a  dense  smoke,  H  can  vary  considerably  even  within  the 
best  spatial  resolution  available  to  the  LCM  (1.5  m).  For  the  program 
of  this  report,  Simpson’s  Rule  was  used  since  It  Is  simple  to  apply  and 
the  theoretical  error  term  la  of  the  fourth  order  in  the  digitization 
rate.  This  Is  the  most  accurate  method  for  numerical  Integration  under 
conditions  of  a  uniform  digitization  rate  and  at  the  same  time  gives  an 
Integrated  value  (In  combination  with  the  trapezoidal  rule)  at  every 
digitized  point.  Table  Cl  lists  the  errors  associated  with  the  numeri¬ 
cal  integration  of  Fig.  Cl  for  ranges  of  offset  from  0-0.9,  signal 
heights  of  1  to  1000  and  an  Infinite  digitization  word  size.  It  is 
evident  that  large  offsets  can  be  tolerated  only  If  the  signal  remains 
relatively  small. 

Cl .2  Digitization  Errors 

The  digitization  errors  have  two  sources  other  than  that  which 
Is  Included  in  the  maserlcal  Integration  error.  These  are  the  finite 
word  length  and  round-off  error.  In  the  calculation  done  here,  the  two 
are  treated  together.  Tables  CII  and  CIII  are  the  errors  obtained  by 
varying  the  word  size,  signal  height  H,  and  the  offset  D.  It  can  be 
seen  that  for  a  certain  combination  of  H  and  D,  there  Is  a  small  range 
of  word  sizes  in  which  the  error  increases  rapidly. 
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C2.0  DISCUSSION  OF  ERRORS 

The  LCM  digitizes  lider  data  with  an  8-bit  word  size.  In  addi¬ 
tion,  the  maximum  observed  relative  change  in  signal  strength  due  to 
red  phosphorus  or  HC  smokes  is  20  times  over  1.5  a.  Assuming,  as  a 
rough  approximation,  that  the  signal  increases  linearly  in  the  inter¬ 
val,  the  maximum  error  due  to  niaerical  integration  occurs  for  an  off¬ 
set  of  0.5  when  the  signal  has  Increased  10  times.  Table  Cl  then  gives 
a  value  of  13Z  for  the  error  due  to  numerical  integration.  For  the 
same  case,  the  digitization  error  becomes  15Z.  In  general,  the  rela¬ 
tive  change  in  the  signal  is  considerably  less  than  20  times  and  errors 
less  than  5-10Z  seem  to  be  more  likely.  The  integration  over  the  total 
signal  will  have  less  error  due  to  less  signal  variation  and  seme  error 
cancellation. 

It  is  clear  from  the  above-mentioned  information  that  computa¬ 
tional  errors  depend  strongly  on  the  signal  variation.  Therefore  the 
Inversion  method  presented  in  this  report  incurs  less  error  than  the 
other  methods  since  transmission  can  be  derived  from  one  Integral  (the 
normalized  Integrated  backscatter)  Instead  of  two  (l.e.  the  normalized 
Integrated  backscatter  and  the  integral  of  extinction  values). 


AGILE 

(AUTOMATICALLY  GENERATED  INVERSION  OF  THE  LIDAR  EQUATION) 
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•••uiptloa.  contained  In  tb.  alngla-acattar  Ildar  aquation  and  In  tha  aaataptlona  contalnad  In  tha  alnpla-acattar  Ildar  aquation  and  In  the 

invaralan  net hod  are  outlined  aa  la  their  relative  Importance.  Special  lnveralon  method  are  outlined  ae  la  their  relative  Importance.  Special 

attention  la  given  to  the  backecatter  veraua  total  extinction  relation,  attention  la  given  to  tha  backecatter  varaua  total  extinction  relation, 

multi-acat taring  and  non-ideal  detector  reeponee.  multl-acatterlng  and  non-ldaal  datactor  reaponaa. 
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